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ABSTRACT 

We investigate the ionization state of the damped Lya system at z = 2.62 toward GB1759+7539 
through an analysis of ionic ratios sensitive to photoionization: Ar°/S + , Fe ++ /Fe + , N + /N°, A1 ++ /A1 + . 
Approximately half of the metals arise in a mostly neutral velocity component with HI/H > 0.9, based 
on Fe ++ /Fe + < 0.013. In contrast, the remaining half exhibits Fe ++ /Fe + « 0.3 indicative of a partially 
ionized medium with HI/H w 0.5. These conclusions are supported by the observed N + /N°, A1 ++ /A1 + , 
and Ar°/Si + ratios. 

We assess ionization corrections for the observed column densities through photoionization models 
derived from the CLOUDY software package. In the neutral gas, the ionization corrections are negligible 
except for Ar°. However for the partially ionized gas, element abundance ratios differ from the ionic 
ratios by 0.1 - 0.3 dex for (Si+, S+, Ni+, Al+)/Fe+ ratios and more for (N°, Ar°)/Fe+. Independent 
of the shape of the photoionizing spectrum and assumptions on the number of ionization phases, these 
ionization corrections have minimal impact (< 0.1 dex) on the total metallicity inferred for this damped 
Lya system. Measurements on the relative elemental abundances of the partially ionized gas, however, 
have a greater than «0.15 dex uncertainty which hides the effects of nucleosynthesis and dust depletion. 

We caution the reader that this damped system is unusual for a number of reasons (e.g. a very low 
Ar°/S + ratio) and we believe its ionization properties are special but not unique. Nevertheless, it clearly 
shows the value of examining photoionization diagnostics like Fe ++ /Fe + in a larger sample of systems. 

Subject headings: galaxies: abundances — galaxies: chemical evolution — quasars : absorption lines — 
nucleosynthesis 



1. INTRODUCTION 

Through echelle optical spectroscopy of the damped Lya 
systems - quasar absorption line systems with A (HI) > 
2 x 10 20 cm -2 - one studies the chemical abundances of 
high redshift galaxies in a manner analogous with UV spec- 
troscopy of the Galactic ISM. These observations yield 
ionic column density measurements of a series low-ions 4 
including Si+, Fe+, Ni+, Zn+, S+, and N° (e.g. Lu et 
al. 1996; Prochaska & Wolfe 1999; Molaro et al. 2000; 
Prochaska et al. 2001), often at a greater accuracy than 
measurements of the Galactic ISM. By studying these chem- 
ical abundances, one gains insight into the nucleosynthetic 

1 http:/ /kingpin. ucsd.edu/~hiresdla 

2 Visiting Astronomer, W.M. Keck Telescope. The Keck Observa- 
tory is a joint facility of the University of California and the Califor- 
nia Institute of Technology. 

3 Hubble fellow 

4 The term low-ion refers to the dominant ionic state of a given 
element in an HI region. 



enrichment, dust properties, and metallicity of the damped 
Lya systems and thereby the high redshift protogalactic 
population (e.g. Prochaska & Wolfe 2002). 

To date, most damped Lya abundance studies have pro- 
ceeded under the assumption that the gas has a very low 
ionization fraction, 



In this case, the low-ion species represent the dominant 
ionization states and ionization corrections for elemen- 
tal abundances are small, e.g. Y(Si) ~ Y(Si + ). There 
is a good theoretical basis for adopting this assumption. 
Prochaska & Wolfe (1996) performed a series of radia- 
tive transfer calculations and predicted x < 10% for sys- 
tems with A(HI) > 10 20 cm~ 2 . Furthermore, the authors 
compared observed A1 ++ /A1 + and Si 3+ /Si + ratios for the 
damped Lya systems against a series of photoionization 
models derived from the CLOUDY software package (v. 
95; Ferland 2001) and concluded it has an ionization frac- 
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tion x < 50%. Viegas (1995) reached similar conclu- 
sions for systems with iV(HI) > 10 21 cm~ 2 but warned 
that ionization corrections could be important in systems 
with iV(HI) w 10 20 cm" 2 . More recently, Howk & Sem- 
bach (1999) argued that photoionization could imply sig- 
nificant corrections for some ionic species in the damped 
Lya systems, particularly if there were a significant local 
source of ionization (e.g. recent star formation). Vladilo 
et al. (2001) have extended this analysis further and con- 
cluded that ionization corrections are a concern for only 
a few of the elements (e.g. Al) observed in the damped 
systems. Finally, Izotov, Schaerer, & Charbonnel (2001) 
have presented a specific model for the damped Lya sys- 
tems which suggests large ionization corrections for a num- 
ber of elements. This model predicts O /Si + ratios, how- 
ever, which are much lower than those typically observed 
in the damped Lya systems. 

All of these studies have presented theoretical arguments 
which require idealized physical conditions (e.g. constant 
volume density, plane-parallel geometry), unknown quan- 
tities (e.g. shape, flux of the ionizing radiation), and un- 
certain physics (e.g. dielectric recombination rates). To 
overcome many of these uncertainties, one can empirically 
determine the ionization state of the gas by measuring 
adjacent ions of several elements (e.g. Fe ++ /Fe + ). Ob- 
servationally, this requires an analysis of transitions like 
Si III A1206, Fe III A1122, and Nil A1083 all of which 
arise in the Lya forest and typically at observed wave- 
length A < 4000A. Unfortunately, we have not tended to 
acquire these observations because of the limited wave- 
length coverage of the HIRES spectrograph in one setting 
and its poor sensitivity below 4000 A. Therefore, the large 
damped Lya databases provided by the Keck community 
(Lu et al. 1996; Prochaska & Wolfe 1999; Pettini et al. 
1999; Prochaska et al. 2001) have not directly addressed 
ionization corrections for the damped Lya abundances. 
Similarly, other studies in the literature (e.g. Fan & Tytler 
1994) have only focused on ratios like C 3+ /C + which pro- 
vide a less precise evaluation of photoionization. 

In this paper, we examine the photoionization of a sin- 
gle damped Lya system through observations of a series of 
ions sensitive to the ionization state of the gas. We present 
a case study of the z a b s = 2.62 damped Lya system toward 
GB1759+7539 (z em = 3.05, V«17). This sightline has 
been previously analysed by Prochaska & Wolfe (1999), 
Outram, Chaffee, & Carswell (1999), and Prochaska et al. 
(2001), but none of these studies focused on the ionization 
state of the damped Lya system. Here, we present obser- 
vations of a series of metal-line transitions which indicate 
this damped Lya system is partially ionized. Before pro- 
ceeding, we emphasize that this damped Lya system was 
selected on the basis that it exhibits several transitions 
suggestive of photoionization (e.g. strong Si III 1206 ab- 
sorption) . Therefore, the system may not be representative 
of the damped Lya population. In fact, there are at least 
four peculiar characteristics of this damped system: 

1. It is the only known damped system with a 
significantly sub-solar Ar°/Si + ratio (Prochaska & 
Wolfe 2002). 



3. It is one of the few systems where the Al III profile 
shows significant departures from the low-ion 
profiles (Wolfe & Prochaska 2000). 

4. It is one of the few systems where the C II* profile 
does not closely trace the low-ion profiles (Wolfe, 
Prochaska, & Gawiser 2002). 

Interestingly, these characteristics are in much better gen- 
eral agreement with gas probed in the Galactic ISM than 
with other damped systems (e.g. Savage and Scmbach 1996). 
Irrespective of the peculiar nature of the damped system 
toward GB1759+7539, our analysis emphasizes the impor- 
tance of observations of Si III, Fe III, N II, Ar I and other 
transitions for assessing the impact of photoionization on 
the chemical abundances one derives from low-ion transi- 
tions. 



Table 1 
ATOMIC DATA 



Transition 


A 




/ 


Rcf 


SIII 1012 


1012 


5020 





03550000 


1 


Arl 1048 


1048 


2199 





26280000 


3 


Arl 1066 


1066 


6600 





06747000 


3 


Nil 1083 


1083 


9900 





10310000 


1 


Felll 1122 


1122 


5260 





16200000 


3 


NI 1134a 


1134 


1653 





01342000 


1 


NI 1134b 


1134 


4149 





02683000 


1 


NI 1134c 


1134 


9803 





04023000 


1 


Fell 1143 


1143 


2260 





01770000 


2 


Sill 1193 


1193 


2897 





49910000 


1 


NI 1199 


1199 


5496 





13280000 


1 


NI 1200a 


1200 


2233 





08849000 


1 


NI 1200b 


1200 


7098 





04423000 


1 


Silll 1206 


1206 


5000 


1 


66000000 


1 


SII 1250 


1250 


5840 





00545300 


1 


CII 1334 


1334 


5323 





12780000 


1 


CII* 1335 


1335 


7077 





11490000 


1 


Nill 1370 


1370 


1310 





07690000 


5 


SilV 1393 


1393 


7550 





52800000 


1 


SilV 1402 


1402 


7700 





26200000 


1 


Sill 1526 


1526 


7066 





12700000 


9 


CIV 1548 


1548 


1950 





19080000 


1 


CIV 1550 


1550 


7700 





09522000 


1 


Fell 1608 


1608 


4511 





05800000 


17 


Fell 1611 


1611 


2005 





00136000 


16 


A1II 1670 


1670 


7874 


1 


88000000 


1 


Nill 1709 


1709 


6042 





03240000 


7 


Nill 1741 


1741 


5531 





04270000 


7 


Sill 1808 


1808 


0130 





00218000 


11 


A1III 1854 


1854 


7164 





53900000 


1 


Znll 2026 


2026 


1360 





48900000 


13 


CrII 2066 


2066 


1610 





05150000 


13 



References. — Key to References — 1: Morton (1991); 
2: Howk et al. (2000); 3: Morton (2002); 

5: Fedchak & Lawler (1999); 7: Fedchak, Wiese, & Lawler (2000); 
9: Schectman, Povolny, & Curtis (1998); 10: Bergeson et al. (1996); 
11: Bergeson & Lawler (1993); 13: Bergeson & Lawler (1993b); 
17: Bergeson, Mullman, k, Lawler (1996b) 

Note. — See http://kingpin.ucsd.edu/~hircsdla for the most 
current and complete version of this table. 



2. It is one of the few systems where the low- ions do 
not uniformly track one another in velocity space. 
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Fig. 1. — Velocity profiles of low-ion transitions from the z a \y a = 2.62 damped Lyoj system toward GB1759+7539. In the figure, v = km/s 
corresponds to z = 2.62561. Most of the transitions are overplotted with a profile fit derived with the VPFIT software package. The marks 
(arrows and x's) at the top of each panel indicate all of the velocity components included in the fits. The arrows designate the components 
included in the analysis of a specific ion whereas the x's have been excluded. Known blends have been dotted out in the figure. 



2. IONIC COLUMN DENSITIES 

The quasar GB1759+7539 has been observed by at least 
three groups (Prochaska & Wolfe 1999; Outram, Chaf- 
fee, & Carswell 1999; Prochaska et al. 2001) with the 
HIRES spectrograph (Vogt et al. 1994) on the Keck I 
telescope. We have combined these data to produce a spec- 
trum with nearly continuous wavelength coverage from 
A = 3400 - 7590A. The Prochaska & Wolfe (1999) and 
Prochaska et al. (2001) data sets were combined with- 
out coadding because the spectral coverage did not over- 
lap. Therefore, the resolution of the original observa- 
tions was preserved at 6 km/s and 8 km/s respectively. 
Finally, we included the data Outram, Chaffee, & Car- 
swell (1999) to fill in the inter-order gaps longward of 
« 5000 A. In terms of the damped Lya system at z = 2.625 
(logiV(HI) = 20.761 ± 0.007; Outram et al. 1999), the 



observations present measurements of transitions ranging 
from S III 1012 to Cr II 2066. 

In the following sub-sections we derive ionic column 
densities for the majority of transitions identified for this 
damped system. We have performed a line-profile fitting 
analysis using the VPFIT software package 5 kindly pro- 
vided by R. Carswell and J. Webb. This least-squares rou- 
tine minimizes a \ 2 matrix of multiple component Voigt 
profile fits to observed line profiles. Each velocity compo- 
nent is described by three parameters - column density, 
redshift, and Doppler parameter (b- value) - which can se- 
lectively be 'tied' together from ion to ion. It should be 
noted that line blending can lead to significant uncertain- 
ties in the Doppler parameter and sometimes redshift of 
various components. On the other hand, our experience 

5 htt p : / / www . ast .cam . ac . uk/ ~ r f c/ vpfft . ht ml 
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is that the column densities remain reasonably robust for 
unsaturated profiles. It is important to note that we in- 
clude solutions for a number of heavily saturated profiles. 
These components are easily identified by their very large 
reported uncertainties and they have not been included in 
the analysis. Furthermore, some velocity components have 
large uncertainties because of a degeneracy in the solution 
(e.g. two velocity components which nearly overlap). We 
have taken care to insure that these line blending effects 
have not influenced our photoionization analysis. Finally, 
Table 1 presents the atomic data used for the metaldine 
transitions analysed in this paper. 

2.1. Low-Ion Profiles 

Figure 1 presents velocity profiles for the most impor- 
tant low-ion transitions in our analysis. The velocity pro- 
files are comprised of two main features at v m and 
—25 km/s where v = km/s corresponds to the redshift 



z = 2.62561. Contrary to the majority of damped Lya 
systems, it is obvious from the figure that the relative 
abundances of these two components are not constant from 
ion to ion. The differences are less severe than the visual 
impression, but are greater than 0.1 dex (i.e. > 3cr) in 
several cases. This is a rare occurrence in the damped 
Lya systems (see Prochaska et al. 2001) and it raises the 
likelihood that this system is peculiar. 

Table 2 presents the full VPFIT solution and errors for 
the low-ion profiles derived with the VPFIT software pack- 
age. In Figure 1, most of the transitions are overplotted 
with the profile fit and the marks (arrows and x's) at the 
top of each panel indicate all of the velocity components 
included in the fits. The arrows designate the components 
included in the analysis of a specific ion whereas the x's 
identify excluded components. In general, all of the low- 
ions were fit with 3-7 components with tied redshifts. In 
contrast to our previous experience with VPFIT solutions 
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Table 2 

IONIC COLUMN DENSITIES: Low-Ions 



ID 


%abs 


a(z) 


Vref 


Ion 


\ogN 


a{N) 


b 


a(b) 






(io- 5 ) 


(km/s) 




(cm" 2 ) 




(km/s) 




1 


2.624329 


4.1 


-105 


C+ 


13.747 


0.514 


20.41 


12.68 










A1+ 


12.677 


0.093 


29.23 


5.28 


2 


2.624360 


0.3 


-103 


Si+ 


13.452 


0.010 


10.67 


0.32 










C+ 


14.160 


0.237 


11.24 


2.78 










A1+ 


12.428 


0.164 


11.48 


1.65 


3 


2.624610 


0.4 


-82 


Si+ 


13.164 


0.042 


5.84 


0.65 










C+ 


14.042 


0.082 


6.37 


1.68 










A1+ 


12.008 


0.093 


4.94 


0.77 


4 


2.624832 


1.3 


-64 


Si+ 


13.505 


0.365 


8.76 


1.55 










Fe+ 


13.453 


0.136 


6.32 


1.04 










C+ 


14.794 


1.180 


7.21 


5.14 










A1+ 


12.383 


0.409 


7.77 


1.82 


5 


2.625035 


0.8 


-47 


Si+ 


14.902 


0.050 


12.40 


1.69 










Fe+ 


14.331 


0.039 


11.03 


1.12 










N° 


14.401 


0.038 


13.72 


1.19 










A1+ 


13.196 


0.119 


14.56 


4.87 










Zn+ 


11.732 


0.153 


7.69 


3.42 










Ni+ 


13.149 


0.040 


15.48 


1.47 










S+ 


14.421 


0.065 


12.72 


2.42 


6 


2.625289 


0.3 


-26 


Si+ 


15.129 


0.024 


8.75 


0.45 










Fe+ 


14.700 


0.016 


11.38 


0.38 










C+ 


16.212 


8.658 


11.10 


23.77 










N° 


14.789 


0.017 


14.08 


0.58 










A1+ 


13.537 


1.752 


5.32 


4.91 










Zn+ 


12.221 


0.060 


11.06 


2.04 










Ni+ 


13.463 


0.019 


12.19 


0.59 










Cr+ 


12.991 


0.094 


13.78 


3.83 










S+ 


14.786 


0.035 


10.78 


0.79 










Ar° 


13.307 


0.028 


9.26 


1.04 


7 


2.625612 


0.1 


+0 


Si+ 


15.057 


0.014 


9.77 


0.37 










C+ 


17.831 


1.661 


8.07 


4.93 










N° 


14.627 


0.018 


5.93 


0.26 










A1+ 


14.132 


0.622 


9.25 


2.24 










S+ 


14.740 


0.126 


7.62 


1.34 










Ar° 


13.432 


0.021 


8.49 


0.67 


8 


2.625644 


0.2 


+3 


Fe+ 


14.473 


0.010 


8.46 


0.21 










Ni+ 


13.208 


0.019 


7.08 


0.51 










Cr+ 


12.826 


0.110 


7.89 


2.73 


9 


2.625914 


0.3 


+25 


Si+ 


14.216 


0.067 


4.29 


0.41 










Fe+ 


13.602 


0.019 


7.80 


0.47 










A1+ 


12.501 


0.157 


3.48 


1.15 










Ni+ 


12.498 


0.097 


9.41 


3.10 










S+ 


14.479 


0.540 


23.68 


32.07 


10 


2.625975 


8.1 


+30 


Si+ 


13.240 


0.568 


8.87 


4.92 










c+ 


14.161 


0.542 


7.22 


4.70 










A1+ 


12.087 


0.562 


7.65 


4.22 


11 


2.626142 


19.9 


+44 


Si+ 


12.784 


0.533 


17.23 


9.25 










c+ 


13.499 


1.107 


14.30 


16.84 










A1+ 


11.560 


0.479 


18.75 


12.47 



"Velocity relative to z=2.62561 

Note. — Portions of the C II 1334 and Al II 1670 transitions arc heavily saturated in some velocity regions and column densities measured 
for these transitions are highly uncertain. 
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Table 3 

IONIC COLUMN DENSITIES: Intermediate-Ions 



ID 


Z u B 
^abs 


do -5 ') 


u rel 

(km /si 


Ion - 


log N 


rr(N\ 


b c 
(km /si 


rr(h) 


1 


9 fi9A^fiO 


0.0 


— 103 


A1++ 


12.264 


0.049 


11.64 


1.07 


2 


9 fi9AK^9 

Zi. UZi'rfcOUZi 


0.0 


—64 


A1++ 


12.128 


0.145 


30.72 


13.42 










N+ 


ro.yuy 


U. 1UZ 


79 


n nn 










re 


ro.yuy 


U. ±o I 


79 


n nn 
u.uu 


Q 
O 




n n 
u.u 


— t i 


A1++ 


12 922 


n n^i 

U.UOl 


1 ^ 87 


1 32 










N+ 


14 047 


n HQS 
u.uyo 


1 ^ 87 


n nn 
u.uu 










Fe++ 


1 Q1 Q 


u.uou 


1 ^ 87 


n nn 
u.uu 


4 


2.625244 


0.7 


-30 


A1++ 


13.357 


0.022 


16.50 


0.85 










N+ 


14.764 


0.045 


16.50 


0.00 










Fe++ 


14.066 


0.066 


16.50 


0.00 


5 


2.625612 


0.0 


+ 


A1++ 


12.629 


0.157 


5.27 


0.89 










N+ 


13.838 


0.232 


5.27 


0.00 


6 


2.625697 


6.5 


+7 


A1++ 


12.429 


0.285 


10.79 


3.05 










N+ 


13.931 


0.205 


6.81 


2.62 


7 


2.625914 


0.0 


+ 25 


A1++ 


12.105 


0.058 


9.29 


0.70 










N+ 


13.918 


0.034 


9.29 


0.00 


8 


2.626142 


0.0 


+ 44 


A1++ 


11.710 


0.329 


27.48 


26.17 










N+ 


13.473 


0.385 


27.48 


0.00 



"Components 1-3,5,8,9 were fixed to have identical redshifts as components observed in the low-ion profiles. 
^Velocity relative to z=2.62561 

c Componcnts with identical b values were tied together in the fit. In this case, an error is listed for only one component. 



of damped Lya profiles (e.g. Prochaska & Wolfe 1996), 
we found a significant improvement in \ 2 if the Doppler 
parameters were left untied from ion to ion. Even more 
unusual, however, is an observed offset in the v « km/s 
component between the Fe II, Ni II, and Cr II transitions 
and the other low-ion transitions. We found a significant 
improvement in \ 2 if we fit the former transitions with a 
velocity component offset by +3 km/s. We have experi- 
mented with our VPFIT solution and have been able to 
shift the ccntroids of the Ni II and Fe II components by 
« 1 km/s closer to the other ions by introducing another 
velocity component at v ~ 10 km/s. The resulting fit, 
however, is not a significant improvement in x 2 over the 
present solution. We are concerned with this offset, in par- 
ticular its implications for analyses which are extremely 
sensitive to such effects (e.g. Webb et al. 2001). It would 
have only a small impact on the current analysis, however, 
and we proceed under the assumption that all of the gas 
near v — km/s has the same velocity. 

2.2. Intermediate- Ion Profiles 

We present a series of intermediate-ion profiles in Fig- 
ure 2 which includes transitions of Fc ++ , Si ++ , Al ++ , and 
N + . To facilitate comparison with the low- ion transitions, 
we also plot the Si II 1808 transition. With the exception 
of Al III 1854, the profiles arc located within the Lya for- 
est and may be contaminated by coincident Lya clouds 
(we are confident that the dotted regions in Figure 2 are 
blends). For this reason, one may consider these column 
density measurements as upper limits except in the case 
where the profiles are heavily saturated and a precise es- 
timate is impossible. We have fit these transitions with 
the VPFIT software package except Si III 1206 which we 
found to be too saturated to include in the analysis. If 
this ion has a similar component structure to the other 
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Fig. 2. — Same as Figure 1 but for the intermediate-ion transitions. 
For comparison, we also include the Si II 1808 low-ion profile. 
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intermediate-ions, then the majority of components have 
N(Si++) > 10 14 cm- 2 . We adopt this limit throughout. 

Similar to Figure 1, the marks at the top of each panel 
indicate the velocity centroids of the components included 
in the analysis with the arrows identifying the specific com- 
ponents included for a given ion. As described below, it 
is an analysis of these ions along with the low-ions which 
places the tightest constraints on the ionization state of 
the damped Lya system. Therefore, we have attempted 
to fit the intermediate-ions with the line-profile solution 
obtained from the low-ion transitions, i.e., the compo- 
nents are placed at the same velocities as the low-ion solu- 
tion. We achieved a satisfactory fit for the intermediate- 
ion profiles except at v « —26 km/s and v « +8 km/s. 
At»« -26 km/s, the x 2 of the fit was significantly im- 
proved if the velocity component was offset by —4 km/s 
from the respective low-ion feature. We contend this off- 
set is a physical effect and its implications on our pho- 
toionization modeling are discussed below. Meanwhile, 
we had to introduce a new component at v « 8 km/s to 
account for significant absorption at that velocity in the 
intermediate-ions. Table 3 presents the full solution for 
the intermediate-ion profiles. 
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Fig. 3. — Same as Figure 1 but for the high-ion transitions. For 
comparison, we also include the Al III 1854 intermediate-ion profile. 



2.3. High-Ion Profiles 

In Figure 3, we present the profiles and VPFIT solution 
for the high-ion Si IV and C IV transitions. We have tied 
the redshifts of the velocity components for the two ions 
but have allowed the 6-values to vary. We chose to fit the 
profiles independently of the low and intermediate-ions - 



even though several of the velocity components have sim- 
ilar redshifts - because we believe the high-ion gas is spa- 
tially separated from the low-ion gas (Wolfe & Prochaska 
2000). For comparison, we plot an intermediate-ion Al III 
transition to emphasize its likeness to the Si IV profiles. 
The parameters of the solution are listed in Table 4. 
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Fig. 4. — Velocity profiles of metal-line transitions from an ab- 
sorption system associated with the z abs = 2.62 damped Lya sys- 
tem toward GB1759+7539. This system has very simple kinematics, 
suggesting the multiple ionic species arise in a single phase of gas. 
Most of the transitions are overplottcd with a profile fit derived with 
the VPFIT software package (Table 5). Known blends have been 
dotted out in the figure. 



2.4. v ~ —330 km/s Feature 

Finally, we analysed a single absorption feature at v = 
—330 km/s which exhibits significant Si + , Si ++ , Si 3+ , C + , 
C 3+ , and A1++ absorption (Figure 4). We fit all of the 
profiles with a single component except C 3+ whose asym- 
metric shape indicates two components. The centroids of 
all of the ions are at nearly identical velocity including the 
stronger component of the C IV profiles. The column den- 
sities for these ions are very accurately measured (Table 5) 
and they provide a precise examination of the ionization 
state of this gas. 

2.5. Summary 

Table 6 presents an observational summary of various 
ionic ratios and column densities for the majority of ve- 
locity components identified in the metal-line profiles. The 
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Table 4 

IONIC COLUMN DENSITIES: High-Ions 



ID 


%abs 


a(z) 
(10- 5 ) 


v rel 

(km/s) 


Ion 


log N 
( cm~ 2 ) 


<j(N) 


b 

(km/s) 


a(b) 


1 


2.624064 


6.6 


-127 


C 3 + 


13.082 


0.097 


27.54 


5.58 


2 


2.624293 


0.3 


-108 


C 3+ 


13.335 


0.036 


5.27 


0.50 










Si 3 + 


13.078 


0.043 


7.99 


0.32 


3 


2.624515 


0.4 


-90 


C 3+ 


13.997 


0.020 


8.16 


0.48 










Si 3 + 


13.062 


0.474 


9.74 


2.39 


4 


2.624769 


4.4 


-69 


C 3+ 


13.344 


0.287 


9.51 


5.39 










Si 3 + 


12.671 


2.049 


19.34 


28.56 


5 


2.624940 


4.1 


-55 


C 3+ 


13.426 


1.201 


6.74 


4.60 










Si 3 + 


13.224 


0.302 


13.23 


10.95 


6 


2.625097 


2.7 


-42 


C 3+ 


14.084 


0.462 


11.57 


11.85 










Si 3 + 


13.721 


0.156 


9.49 


1.45 


7 


2.625291 


2.4 


-26 


C 3+ 


13.467 


1.444 


6.61 


8.54 










Si 3 + 


13.419 


0.297 


9.62 


4.71 


8 


2.625417 


11.5 


-15 


C 3+ 


13.374 


0.943 


7.92 


8.63 










Si 3 + 


13.020 


1.713 


13.49 


47.17 


9 


2.625601 


5.8 





C 3+ 


13.390 


0.296 


10.82 


5.37 










Si 3 + 


12.830 


1.236 


9.37 


4.35 


10 


2.625880 


2.9 


+22 


C 3+ 


13.199 


0.117 


9.86 


2.53 










Si 3 + 


12.847 


0.199 


12.26 


5.16 


11 


2.626038 


2.8 


+35 


C 3+ 


12.676 


0.341 


7.13 


3.59 










Si 3 + 


12.508 


0.336 


9.39 


4.50 


12 


2.626206 


1.0 


+49 


C 3+ 


12.603 


0.106 


5.01 


1.38 










Si 3 + 


12.341 


0.137 


5.09 


1.02 


13 


2.626430 


0.5 


+68 


C 3+ 


12.847 


0.047 


16.30 


2.08 










Si 3 + 


12.236 


0.065 


18.97 


3.30 



Table 5 

IONIC COLUMN DENSITIES: v = -330km/s Component 



Ion 


%abs 


a(z) 


Vrel 


log N 


a(N) 


b 


a(b) 






(io- 5 ) 


(km/s) 


( cm~ 2 ) 




(km/s) 




C+ 


2.621580 


0.6 


-332 


12.775 


0.035 


4.36 


0.96 


C 3+ 


2.621484 


12.1 


-340 


13.030 


0.876 


6.62 


5.17 


C 3+ 


2.621580 


3.0 


-332 


13.455 


0.337 


5.45 


1.33 


Si+ 


2.621590 


0.5 


-332 


12.387 


0.031 


4.64 


0.85 


Si++ 


2.621584 


0.3 


-332 


12.526 


0.028 


5.00 


0.44 


Si 3 + 


2.621570 


0.1 


-333 


12.956 


0.016 


3.54 


0.15 


A1++ 


2.621571 


0.8 


-333 


11.537 


0.112 


5.39 


2.39 


A1+ 






-333 


< 10.8 
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Table 6 
SUMMARY TABLE 
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0.77 


-0.01 


-0.30 
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-1.05 




%[AW(Si+)] 


0.1 


0.8 


0.9 


22.5 


37.9 




32.2 


4.6 


%{N TO t(S + )} 
%[N TO t(Fc+)} 








15.3 


35.4 




31.9 


17.5 






2.6 


19.8 


46.4 




27.5 


3.7 


%[N TO T(Fe ++ )} 






28.5 


29.2 


40.9 




< 1.4 




%[N TO t(M++)} 


0.8 


4.2 


3.1 


19.3 


52.5 




9.8 


2.9 


%[Ntot(N+)] 






7.8 


10.7 


55.8 




6.6 


8.0 


%[A^ TO t(N )] 






< 3.7 


18.8 


45.9 




31.6 





velocities listed at the top of the Table designate 10 km/s 
bins where all of the gas within 5 km/s of the centroid 
has been summed for the ionic column densities. All ra- 
tios are logarithmic values with the square bracket nota- 
tion indicating values relative to solar abundances, i.e., 
[X/Y] eee Iog(JV(X)/JV(Y)) - log(X/Y) Q . 

3. ANALYSIS OF THE GAS AT V < -50 km/s 

In this section we consider various diagnostics on the 
ionization state of the gas at v < —50 km/s along the 
sightline. It is important to keep in mind that the ab- 
sorption at any velocity may be the result of several com- 
ponents which simply overlap in velocity space. These 
components might each have a unique set of physical char- 
acteristics (e.g. mctallicity, ionization state) and therefore 
comparisons of different ions at a given velocity may be 
open to several interpretations. 

3.1. Gas at v k, —330 km/s 

As revealed by Figure 4, there is an absorption feature 
at v w —330 km/s which exhibits significant absorption 
for many ions: Si+, Si++, Si 3+ , C 3 +, C+, and A1++. The 
kinematic structure is very simple and we fit the profiles 
with a single component (Table 5). In fact, the feature 
is so simple that we are confident the majority of gas is 
spatially well- localized. The Si 3+ , C 3+ , and Al ++ ions 
have significantly higher column densities than their low- 
ion counterparts indicating the gas is highly ionized. Sig- 
nificant Si 3+ and C 3+ column densities can arise from a 



number of physical processes including collisional ioniza- 
tion and photoionization by high energy photons. Because 
the observed Doppler parameters are small, we place a 2a 
upper limit on the temperature independently from the 
C + and Si 3+ measurements: T < 3 x 10 4 K. In turn, this 
indicates the gas is not in collisional ionization equilibrium 
because T > 4 x 10 4 K is required to produce significant 
amounts of C 3+ (c.f. Sutherland & Dopita 1993). This 
naturally suggests the ionization may be dominated by 
photoionization. 

Although this system is well separated by velocity from 
the damped Lya system, its characteristics are important 
for several reasons. First, it is possible that this gas is rep- 
resentative of the gas giving rise to the C IV and Si IV pro- 
files presented in Figure 3. Therefore, the gas responsible 
for the C IV and Si IV profiles makes a negligible contri- 
bution to the low-ion column densities - and presumably 
AT (HI) - of the damped Lya system. A similar conclusion 
has been inferred from the observed differences in the low 
and high-ion profiles of many other damped Lycn systems 
(Lu ct al. 1996; Wolfe & Prochaska 2000). We also note 
that the Si ++ /Si 3+ and Al++/Si 3+ ratios are sufficiently 
small for the v w —330 km/s component that gas with its 
ionization characteristics cannot reproduce both the high 
and intermediate-ion column densities at v « km/s. If 
we attribute the CIV and Si IV profiles at v > —50 km/s 
in Figure 3 to an ionized gas similar to the v w —330 km/s 
component, then we must introduce at least one additional 
phase of gas for the intermediate-ions independent of the 
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low-ions. Therefore, if the intermediate and low-ions have 
distinct origins, this indicates the presence of three unique 
ionization phases at a similar velocity in this damped Lya 
system. 

The absorption feature at v w —330 km/s is further im- 
portant because it poses a difficult challenge to ionization 
models. Because of its simple kinematics and low temper- 
ature, the gas is very likely in a single component with a 
single set of physical conditions. Equilibrium collisional 
ionization models or more exotic turbulent mixing layers 
(e.g. Shull & Slavin 1994) are inconsistent with the tem- 
perature of the gas. If the gas is in equilibrium, therefore, 
then we presume it has been photoionized. 

With the range of ions and very accurate column densi- 
ties observed, we expected to precisely constrain its ioniza- 
tion state. To this end, we performed a series of photoion- 
ization calculations using the CLOUDY software pack- 
age with both an extragalactic UV background spectrum 
(Haardt & Madau 1996) and a series of Kurucz (1988) 
O-star spectra. We varied the HI column density of the 
plane-parallel 'cloud', assumed a range of metallicities, and 
calculated the ionization state at a range of ionization pa- 
rameters, 



U = 



> 



cn H 



(2 x 1(T 5 ) 



J912/IO 



-21.5 



uh/ cm 3 



(2) 



where tih is the volume density of hydrogen, <j> is the sur- 
face flux of ionizing photons with hv > 1 Ryd, and J912 
is the intensity of the incident radiation at 1 Ryd. The 
relation in Equation 2 approximately holds for ionizing 
radiation fields which steeply decline at energies greater 
than 1 Ryd. 

To our surprise, we cannot identify any single-phase 
(e.g. constant density) photoionization model which can 
match the observed ratios of the Si ions (Si + : Si ++ : Si 3+ = 
1:1.3:4). The principal conflict is that models which 
predict iV(Si+) « iV(Si ++ ) always require 7V(Si ++ ) > 
_/V(Si 3+ ) in stark contradiction with our observations. This 
difference between the predicted and observed ratios is 
generally greater than a factor of 10. It holds indepen- 
dently of the shape or flux of the input spectrum or any 
assumptions on the metallicity or H I column density. We 
also note that even if we allow the temperature to be a free 
parameter, it is impossible to reproduce the Si ionic ratios 
through a combination of collisional and photoionization. 
In conclusion, we contend the gas cannot be described by 
a single ionization phase in equilibrium conditions. While 
this could be a common occurrence in quasar absorption 
line systems, it is very surprising for this feature because 
of its simple kinematics. 

Granted the failure of the single-phase equilibrium mod- 
els, we also considered the prospect that the gas arises in 
two distinct phases with very different ionization states. If 
we combine a very neutral component with a highly ion- 
ized component, we approach the observed ratios of Si and 
C but still fail at the > 95% c.l. Furthermore, the best 
fitting model couples an extremely neutral gas with an ex- 
tremely photoionized gas in a mixture of w 2 : 1, which 
would be unlikely given the extremely simple kinematics 
of the observed metal- line profiles. It would take a con- 
trived set of physical conditions to couple two-phases of 
gas with such extreme ionization characteristics at identi- 



cal velocities and both with T < 3 x 10 4 K. Finally, this 
two-phase model predicts an O I column density in conflict 
with our observed upper limit. To this end, we have failed 
to appropriately model the ionization state of this veloc- 
ity component. Although it has no direct bearing on our 
damped Lya analysis because of its very small metal con- 
tent, we interpret this failure as a stern warning for the 
challenges associated with making ionization corrections 
for highly ionized gas. 

Our photoionization models do not accurately predict 
the ionization structure of the —330 km/s complex; fur- 
thermore, collisional ionization equilibrium models can- 
not provide the required ionization since the upper lim- 
its on the kinetic temperature (T < 3 x 10 4 K) are too 
small to allow collisional production of CIV and Si IV in 
great quantities (Sutherland & Dopita 1993). It is pos- 
sible that the gas in this component is not in ionization 
equilibrium. Non-equilibrium ionization can occur, for ex- 
ample, if the cooling time of shock-heated gas is signif- 
icantly longer than the recombination time (Shapiro & 
Moore 1976; Sutherland & Dopita 1993). Models of time- 
dependent ionization of hot gas (e.g., gas hot enough to 
produce significant amounts of the highly-ionized species 
observed in this cloud) generically predict much larger ion- 
ization fractions of highly-ionized species at the tempera- 
tures appropriate for this component than do equilibrium 
models. Thus another potentially valid model for the —330 
km/s gas is that this component was initially shock- heated 
to high temperatures (T > 10 5 K) and has cooled rela- 
tively quickly, with its ionization state lagging its temper- 
ature. While the ionization fractions of C 3+ and Si 3+ can 
be large even to temperatures as low as 10 4 K in such sce- 
narios (Sutherland & Dopita 1993), we note that models 
of non-equilibrium ionization balance are dependent upon 
the initial conditions of the gas. Thus it is difficult to 
determine if the ionization structure of the —330 km/s ab- 
sorption is consistent with this non-equilibrium scenario 
without unfounded guesses at its initial temperature, den- 
sity, and abundance. 

3.2. Gas at v w —105 and —65 km/s 

The low column density gas with v « —105 and —65 km/s 
exhibits significantly different ionization characteristics than 
the v « —330 km/s component. Although the ionic ratios 
(e.g. Si++/Si+, Fc++/Fc+, and N+/N ) all indicate the 
gas is predominantly ionized (x « 1), the C 3+ and Si 3+ 
ions comprise a much smaller fraction of the elemental col- 
umn densities than the v w —330 km/s component. 

The properties of the gas at v « —65, —105 km/s may 
have important implications for the intermediate-ion gas 
observed at v > —50 km/s. Consider the possibility that 
the low and intermediate- ion profiles at v > —50 km/s 
arise from two ionization phases: (i) a neutral gas phase 
and (ii) a highly ionized gas with properties similar to the 
gas at » ~ —65, —105 km/s. In this case, we can esti- 
mate the contribution of this ionized phase to the total 
Fe + column density at v > — 50 km/s from the following 
treatment. At v w —65 km/s, we measure a Fc ++ /Fe + 
ratio of +0.46 dex. If we assume the same ratio applies at 
v « —26 km/s, then the 7V(Fe ++ ) value at v « —26 km/s 
implies N(Fe + ) = 10 14 cnT 2 which is « 20% of the to- 
tal N(Fe + ) value at that velocity. Therefore, an ionized 
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gas component could comprise roughly 20% of the low-ion 
absorption observed at v w —26 km/s. We develop this 
treatment to much greater depth in the following section. 

4. ANALYSIS OF THE LOW-ION REGIONS AT 

v > —50 km/s 

The three velocity components at v « —47, —26, and 
km/s comprise nearly all of the low-ion gas in this damped 
Lya system. In terms of the system's chemical abun- 
dances, therefore, ionization corrections on the column 
densities derived from the low-ion profiles at these veloci- 
ties are the most important. 

4.1. v rts km/s 

It is most instructive to begin with the v « km/s 
component where, initially, there are conflicting photoion- 
ization diagnostics. On the one hand, the Ar°/S + ratio 
is significantly sub-solar which indicates the majority of 
Ar is ionized (e.g. Sofia & Jenkins 1998): [Ar°/S+] « 
—0.7 dex. Therefore, one might expect photoionization 
is significantly affecting the chemical abundances inferred 
from all of the low- ion transitions. On the other hand, we 
have placed a strict upper limit on the Fe ++ column den- 
sity which limits the level of photoionization of this gas. 
The adopted limit of log A(Fe ++ ) < 12.6 dex results from 
the low optical depth observed at v « km/s in the FeUI 
1122 profile (Figure 2). We stress that the absorption 
present at v « km/s is probably due to an intervening 
Lya cloud and therefore the A(Fe ++ ) value is possibly 
lower than 10 12 cm~ 2 . Even adopting our more conserva- 
tive upper limit, Fe++/Fc+ < -1.9 dex (Tables 2,3), this 
indicates the gas has a small ionization fraction. Finally, 
the moderate strength charge-exchange reactions between 
N and H which roughly couple the column densities of 
N°, N + , H°, and H + , allow one to crudely estimate the 
ionization fraction of hydrogen from the observed N + and 
N° column densities. We find N+/(N++ N°) « 0.15 which 
suggests x w 15%, although N can be modestly overionized 
with respect to hydrogen because of its larger ionization 
cross-section. We contend a similar effect explains the sub- 
solar Ar°/S + ratio observed for this predominantly neutral 
gas. 

To perform a quantitative analysis of photoionization, 
we have performed a series of calculations with the CLOUDY 
software package. We caution that the accuracy of these 
calculations are limited by a number of uncertainties asso- 
ciated with the physical nature of the absorbing gas and 
the incident ionization field. We consider here a series of 
single-phase, constant density models with varying ioniza- 
tion parameter, the Haardt-Madau spectrum, and assume 
AT (HI) = 10 20 cm~ 2 which is roughly consistent with the 
relative strength of the low-ions and the total H I. In Fig- 
ure 5, we plot the predicted logarithmic abundance of low- 
ion Xj relative to Fe + relative to the intrinsic abundance 
of these two elements, 

[X;/Fe+] = log(A(X i )/A(Fe+)) - log(X/Fe) , (3) 

against the Fc ++ /Fc + ratio. We focus on [Xi/Fe + ] in- 
stead of [Xj/H°] because we are more concerned with the 
effects of photoionization on interpretations of relative el- 
emental abundances [X/Y] instead of metallicity [X/H] in 
the damped Lya systems. Departures of [Xi/Fe + ] from 



zero indicate that photoionization corrections are required 
to calculate accurate elemental abundance ratios from the 
low-ion column densities independent of dust depletion. 
Positive (negative) [Xi/Fe + ] values would imply overesti- 
mates (underestimates) of [X/Fe] from the low-ion ratios. 
Finally, the black dashed line plots the ionization fraction, 
logs, on the same logarithmic scale. 

At Fe ++ /Fe + < — 1.9 dex, we note log a; < — 1 indicating 
the hydrogen gas is less than 10% ionized. Therefore, we 
are confident that the gas at v ~ km/s is predominantly 
neutral. In turn, we find the ionization corrections are 
negligible for Si + , Cr+, and Ni + relative to Fc+ and small 
for S + , Al + , and N°. Only Ar° requires a large ionization 
correction which brings the corrected Ar/S ratio into closer 
agreement with the solar abundance. This gas resembles 
the (i Cen A ISM sightlinc analysed by Sofia & Jenkins 
(1998) which shows similarly over-ionized Ar/S in a gas 
with a relatively large neutral hydrogen fraction. The fact 
that the inferred Ar/S ratio remains sub-solar even after 
the ionization corrections implied by Figure 5 may reflect 
differential depletion but more likely points to an error in 
our adopted ionization spectrum or possibly uncertainties 
in the atomic physics of Ar. In either case, the difference 
reflects an error in our photoionization modeling. 

We emphasize that our conclusions on the ionization 
state of this gas are further supported by the low A1 ++ /A1 + 
and Fe ++ /Al ++ ratios. Our results are also consistent 
with the observed Si++/Si+, Si 3 +/Si+, and C 3+ /C+ ra- 
tios. Furthermore, we have also considered a Kurucz stel- 
lar ionizing spectrum with T = 30000 — 50000 K and more 
complicated two-phase scenarios and have found no qual- 
itative differences in these analyses. In conclusion, the 
gas at v ~ km/s has a low ionization fraction (x -C 1) 
and ionization corrections to the observed low-ion column 
densities are nearly negligible except in the case of Ar°. 
Therefore, the chemical abundances for this gas are accu- 
rately represented by the low-ion species. 

4.2. v « -26 and -47 km/s 

Now consider the gas at v « —47 and —26 km/s which 
exhibits significantly larger N+/N , A1++/A1+, and Fc ++ /Fc+ 
ratios than the gas at v « km/s. In the case of N + and 
Fe ++ , we stress that these transitions lie in the Lya for- 
est and may suffer from contamination by coincident Lya 
clouds. In the following, we proceed under the assumption 
that such contamination is small. Making a crude estimate 
of the ionization fraction from the observed N + /N° ratio 
as above, we find x « 50% assuming N + /N° « —0.1 dex. 
We also note that the Al ++ /Si + ratios are high relative to 
the majority of damped Lya systems (Vladilo et al. 2001). 
Altogether, we expect this gas to be more highly ionized 
than the v ~ km/s component. 

The greatest challenge in accurately assessing the im- 
pact of photoionization in the damped Lya systems is our 
ignorance of the physical nature of the gas 'clouds' re- 
sponsible for the absorption. We have no a priori knowl- 
edge of the metallicity, HI column density, volume density, 
or the underlying nuclcosynthetic or dust depletion pat- 
terns. Furthermore, the absorption at any velocity could 
very well be the superposition of many clouds each with 
a unique set of physical characteristics but coincident ve- 
locity along the sightlinc. As such, there may be several 
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Fig. 5. — Predicted logarithmic abundance of low-ion Xj relative to Fe+ relative to the intrinsic abundance of these two elements (Equation 3) 
against the Fe ++ /Fe + ratio. Departures of [Xi/Fe + ] from zero indicate that photoionization corrections are required to calculate accurate 
elemental abundance ratios from the low-ion column densities independent of dust depiction. Positive (negative) [Xj/Fe + ] values would imply 
overestimates (underestimates) of [X/Fc] from the low-ion ratios. The black dashed line plots the ionization fraction, logrr, on the same 
logarithmic scale as [Xj/Fe"*"]. For this analysis, we have assumed the Haardt-Madau spectrum, 7V(HI) = 10 20 cm -2 , and [Fc/H] = — 1. All 
of these quantities were derived with the CLOUDY software package. 



clouds with different ionization states giving rise to the 
various metal-line profiles. Our approach is to consider 
simplest models first and then qualitatively assess the ef- 
fects of more sophisticated scenarios. 

We consider the two following cases: (1) the gas arises 
in a single cloud with one ionization phase; and (2) the 
gas arises in two gas clouds with very different degrees 
of ionization and possibly metallicity. Before proceeding, 
we present a line of circumstantial evidence that the ion- 
ization corrections to the elemental abundances might be 
small for this gas. In Table 6, we note that the relative 
abundances of the low-ion species of S, N, Ar, Si, Fe, 
Ni, and Cr are relatively uniform over all of these low- 
ion components. In all cases, the ratios vary by less than 
0.2 dex from component to component. If one adopts the 
unfounded assumption that the gas in the entire damped 
Lya system shares similar elemental abundances, then one 
would infer that the ionization corrections are small. Con- 
versely, significant ionization corrections might require large 
abundance variations, possibly from different chemical en- 
richment histories or dust depletion patterns. 

4.2.1. Single Phase Model 



The simplest scenario is that the gas is contained within 
a single, constant density cloud with one ionization phase 6 . 
If the cloud is subjected to an external radiation field, 
there will be a radial variation in the degree of ionization 
of the gas; the outer layers are highly ionized and the in- 
ner regions are primarily neutral. With the exception of 
Al ++ (e.g. Prochaska & Wolfe 1996), the intermediate- 
ions arise in the outer layers of the cloud because few high 
energy photons reach the inner regions. In this important 
respect, the single-phase model mimics the two-phase sce- 
narios presented next. Because the N + /N°, Fe ++ /Fe + , 
and Fe ++ /Al ++ ratios are similar for the gas at v « — 26 
and —47 km/s, we consider the two components together. 
The uncertainties in our photoionization modeling pre- 
sumably exceed any differences in the ionic column densi- 
ties of the two velocity components. 

Our analysis relies on simulations drawn from the CLOUDY 
software package. Figure 6 presents the predicted column 
densities N(X.) for a series of ionization parameters as- 
suming N(HT) = 10 20 cm -2 , a metallicity [Fe/H] = -1, 
and the Haardt-Madau EUVB spectrum at z = 2.5. The 
HI column density and metallicity are reasonably consis- 

6 This treatment is identical to the scenario envisioned by Prochaska 
& Wolfe (1996) and the model designated HI by Vladilo et al. 
(2001). 
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Fig. 6. — Predicted column densities N(X) for a series of ionization parameters assuming Af(HI) = 10 20 cm -2 , a mctallicity [Fc/H] = — 1, 
and the Haardt-Madau EUVB spectrum at z = 2.5. The vertical dashed lines indicate the observed ratios of Fe ++ /Fe+, N+/N , and Ar°/S+ 
and we note that the A1++/A1+ upper limit is consistent requires logt/ < —2 dex. We find that all four of these ionization diagnostics are 
consistent with a single ionization phase and logC/ ~ — 3 dex. 



tent with the total quantities inferred from the entire sys- 
tem. We emphasize that at AT (HI) > 10 20 cm -2 , the exact 
iV(HI) value that one assumes has little impact on the 
single-phase scenarios. We also note that the N(X) values 
presented in Figure 6 are less important than the relative 
values of the pairs of lines with identical color and symbols 
(e.g. Fe ++ /Fe + ). The vertical dashed lines indicate the 
observed ratios of Fe++/Fe+, N+/N , and Ar°/S+ and the 
A1++/A1+ upper limit requires logt/ < — 2 dex. We find 
that all four of these ionization diagnostics are consistent 
with a single ionization phase where log U « — 3 dex. The 
agreement is not perfect, but it is acceptable given the 
number of assumptions and uncertainties in the photoion- 
ization modeling. For example, a slightly softer ionizing 
spectrum would help bring the predicted Ar°/S + ratio into 
better agreement with the observations at log U « — 3 dex. 

The implications of the single phase model with log U « 
—3 are as follows. First, we derive an ionization frac- 
tion x w 50% which indicates the gas is partially ionized. 
Second, if the metallicity of the velocity components at 
v > —50 km/s is constant, then the v w —26 km/s com- 
ponent may exhibit a larger Af(H) value than the gas at 
v w km/s. Third, the degree of photoionization is signif- 
icant enough to imply ionization corrections to the relative 
abundances derived solely from the low-ion species. Re- 
ferring to Figure 5, one notes that as Fc ++ /Fe+ — > dex 



several of the low-ion ratios [Xi/Fc + ] require modest ion- 
ization corrections. With the exception of Ar°, the cor- 
rections are generally small (< 0.3 dex) but significant 
compared to the typical statistical errors in these mea- 
surements. The vertical black dash-dot line in Figure 5 
indicates the approximate observed Fe ++ /Fe + value of 
—0.5 dex. For this value, the low-ion ratios overestimate 
Si and Ni relative to Fe by ~ 0.1 dex while Al, S, and N are 
underestimated by « 0.2 dex. These ionization corrections 
are relatively small, but can have important implications 
for the relative elemental abundances of the damped Lya 
systems. For example, correcting the S and Si abundances 
gives [S/Si] and [Si/Fe] values for the u w -26 km/s com- 
ponent which differ by w 0.3 dex from the v w km/s 
gas. One concludes, therefore, that the relative elemental 
abundances may differ significantly within this damped 
Lya system. 

We have also considered a single-phase model with the 
Kurucz stellar ionizing spectrum and a range of stellar ef- 
fective temperatures. These models imply similar correc- 
tions; compare against Figure 7 which plots aT = 40000 K 
model with AT (HI) = 10 20 cm~ 2 . The principle differences 
are that [Cr + /Fe + ] shows a large offset from solar and that 
[S + /Fe + ] is positive at high Fe ++ /Fc + values and signifi- 
cantly negative at Fe ++ /Fe + < — 2 dex. These differences 
arise from the softer Kurucz ionizing spectra (e.g. a signif- 
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Fig. 7. — Predicted logarithmic abundance of low-ion Xj relative to Fe + relative to the intrinsic abundance of these two elements against 
the Fe ++ /Fe + ratio. The analysis and presentation is identical to Figure 5 except the ionizing spectrum is a Kurucz T = 40000 K stellar 
model. 



icant fraction of S is in S° at low ionization potential) and 
they highlight the uncertainties in our models related to 
the assumed input spectrum. In theory, one might observe 
enough ions to discriminate between various ionizing spec- 
tra but in practice one must allow for these uncertainties. 

4.2.2. Two-Phase Model 

Although we found that the single-phase scenario gave 
reasonable agreement to our set of photoionization diag- 
nostics, we also wish to investigate the effects of a two- 
phase scenario. In particular, we want to determine if the 
two approaches yield different estimates for the ionization 
corrections (e.g Vladilo et al. 2001). We will show that the 
two models predict similar corrections but with important 
differences for S + and Ar°. 

Consider the following scenario where the low and 
intermediate-ions arise in two parcels of gas with unique 
ionization characteristics. As discussed in § 3.2, it is pos- 
sible that gas with properties similar to the ionized gas at 
v w —65 km/s is entirely responsible for the intermediate- 
ions but also contributes a significant percentage of the 
low-ion column densities. In fact, we estimated an ion- 
ized component could contribute « 20% of the Fe + gas 
observed at v « —26 km/s. This scenario resembles the 
physical conditions along the sightlinc to HD 93521 in 
the Galactic ISM, where Sembach et al. (2000) describe 
the impact of the Warm Ionized Medium on measuring 



elemental abundances. For the damped system, a multi- 
phase interpretation is further supported by an observed 
offset between the low and intermediate-ion gas at v « 
—26 km/s. Examining Tables 2 and 3, one notes that 
the component which best describes the low-ion gas at 
v w —26 km/s is offset by w 4 km/s from the centroids 
of the Al ++ , N + , and Fe ++ component. The difference 
is relatively small and may be the result of a system- 
atic error, but we consider it more likely that the gas at 
v ss — 26 km/s is best described by at least two compo- 
nents with very different ionization fractions 7 . In any case, 
because the offset is small one must consider the impact of 
the ionized gas on the observed low-ion column densities 
and we treat the two clouds as if they had identical veloc- 
ities. This treatment is similar to the scenario introduced 
by Howk & Sembach (1999) and extended by Vladilo et 
al. (2001). 

Consider the following model which assumes the gas 
at v f=3 —26 km/s is comprised of two components: (i) a 
very neutral component (e.g. Fe ++ /Fe + <?C 1) and; (ii) an 
ionized component with identical relative gas-phase abun- 
dances that contributes all of the intermediate-ions ob- 
served at v ?3 —30 km/s. Even though the latter com- 
ponent is significantly ionized, it may still exhibit signifi- 

7 Unfortunately, the gas at v ~ —47 km/s is not precisely localized 
in velocity space and can be well described by a single or multiple 
clouds. 
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Fig. 8. — Offsets from the intrinsic abundances of low-ion Xj relative to Fe+, [X;/Fe+], against the fraction of Fe+ arising in the neutral 
component /jv(Fe + )of our two-phase models. Similar to the single- phase models, the implied ionization corrections are small for the majority 
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cant column densities of several low- ion species, (e.g. S + , 
Si + , and Al + ) and, therefore, impact the relative elemen- 
tal abundances and metallicity derived from the low-ion 
transitions. 

To illustrate the ionization corrections for this two-phase 
model we introduce a new parameter, 
f N (Xi) = (X l ) neutra i/(X i ) tot ai, which is the fraction of 
low-ion Xi that is attributed to the neutral gas component. 
We have used CLOUDY calculations similar to those pre- 
sented in the previous section to calculate the ionization 
state of the ionized component but assuming the HI col- 
umn densities are lower than 10 20 cm~ 2 . In this case the 
ionized component contributes a small fraction of the to- 
tal HI column density. We constrain the two-phase model 
by requiring that it yields a Fe ++ /Fc + ratio of —0.5 dex. 
Because this constraint refers to a relative abundance in- 
stead of the absolute Fe ++ and Fe + column densities, our 
analysis is not very sensitive to the A(HI) value or metal- 
licity assumed for the ionized component. Changing the 
AT(HI) value for the ionized component by 0.5 dex can be 
important, however, and we consider the effects below. 

We first present results assuming iV(HI) = 10 19 cm~ 2 in 
the ionized component and adopting the Haardt-Madau 
spectrum. Figure 8 plots the offset from the intrinsic abun- 
dances of low-ion Xj relative to Fe + , [Xi/Fe + ], against the 
fraction of Fe + arising in the neutral component /at(Fc + ). 
As /iv(Fe + ) increases in the figure, the ionized component 
contributes less Fe + and therefore there is a compensatory 
increase in the ionization state of this gas to match the 
observed Fe ++ /Fe + ratio. At low /Ar(Fe + ) values, one ap- 
proaches the single-phase model presented in the previous 
sub-section. Note that one derives a qualitatively similar 
set of results by examining /jv(N°) and constraining the 
predicted N+/N ratio by the observed value. 

Similar to the single-phase model, the implied ioniza- 
tion corrections are small for the majority of elements al- 
though there are several quantitative differences. First, 
the offsets for [S + /Fe + ] are generally small or slightly 
positive whereas the single-phase model predicted neg- 
ative offsets approaching —0.2 dex. The difference can 
be attributed to the difference in optical depth of the 
single-phase model and the ionized component in the two- 
phase scenario where we have assumed iV(HI) = 10 20 cm~ 2 
and 10 19 cm -2 respectively. In the former case, the self- 
shielding from H I is large enough that S is much less ion- 
ized relative to Fe. Second, the offsets for Ar°/Fe + are 
small at all but the lowest /jv(Fe + ) values. This is be- 
cause we have assumed Ar is entirely in atomic form in 
the neutral component. If this two phase model is an ac- 
curate representation of the gas, the neutral component 
must have a sub-solar Ar°/S + ratio similar to the the 
neutral gas at v « Okm/s. Finally, we note large cor- 
rections for A1+, Si+, and S+ at log /at(Fc+) > -0.07 
which corresponds to a very large ionization parameter 
for the ionized component. At very high levels of ioniza- 
tion, the Si + , S + , and Al + ions exhibit significantly greater 
abundances relative to Fe + . We doubt, however, that the 
gas considered here has log /jv(Fe + ) > —0.07. As noted 
above, if the absorption at v w —65 km/s is representa- 
tive of the ionization component in our two-phase model, 
then log /w(Fc + ) < — 0.1 . Furthermore, it is very unlikely 
(if not nonphysical) to assume that the neutral compo- 



nent contributes no Fc ++ gas. Allowing some ionization 
of Fe + to Fe ++ in the neutral component drives the effec- 
tive /iv(Fe + ) value even lower. 

We have also performed a two-phase analysis assum- 
ing a series of Kurucz stellar spectra. Figure 9 presents 
the results assuming AT (HI) = 10 19 cm- 2 , T = 40000 K, 
[Fe/H] = -I, and Fc++/Fc+ = -0.5 dex. The results 
are qualitatively similar to those from the Haardt-Madau 
spectra (Figure 8): Si+, A1+, S+, Cr+, and Ni+ all re- 
quire small negative ionization corrections relative to Fe + 
while N° and Ar° require large positive corrections. In 
contrast with the Haardt-Madau results, all of the correc- 
tions converge to dex as /Ar(Fe + ) — > I. For a stellar 
atmosphere with T — 40000 K, one predicts significantly 
larger Fe ++ /Fe + ratios than the Haardt-Madau spectrum, 
primarily due to the softer spectral slope of the stellar 
flux. In fact, adopting a higher effective temperature (e.g. 
50000 K) for the stellar atmosphere brings the results into 
better qualitative agreement with the Haardt-Madau mod- 
els. Finally, we have also considered different values of 
AT (HI) in the ionized component than 10 19 cm~ 2 . In gen- 
eral, lower AT(HI) values lead to slightly larger [Xi/Fe + ] 
values and higher A(HI) values imply smaller corrections. 
The differences are typically smaller than 0.1 dex, how- 
ever, and are therefore on the level of the uncertainties 
inherent to our photoionization modeling. 

To summarize, the two-phase scenario implies small but 
important corrections to the low-ion ratios. We favor this 
model over the single phase model because of the like- 
lihood that the intermediate and low-ion components at 
v w —26 km/s are not exactly aligned. Furthermore, we 
have identified ionized gas near the damped Lya system 
(at v « —65 km/s) which exhibit properties consistent 
with the ionized component of our two-phase model. In- 
terestingly, the ionization corrections of the two-phase sce- 
nario lead to relative abundances which are less conflict- 
ing with the v « km/s gas. The challenge of the two 
phase scenario in comparison with the single phase model, 
however, is that it introduces an additional gas compo- 
nent with another set of uncertainties related to HI col- 
umn density, metallicity, geometry, etc. These uncertain- 
ties prohibit one from determining the //v(Fe + ) value and 
thereby precise ionization corrections. 

5. DISCUSSION 

In the previous sections, we carefully considered the ion- 
ization state of the gas comprising the damped Lya sys- 
tem at z = 2.62 toward GB1759+7539. Qualitatively, we 
demonstrated that the gas in one major velocity compo- 
nent is primarily neutral (x < 0.1) while the remaining 
gas is partially ionized (x ~ 0.5). In both cases, signifi- 
cant ionization corrections for Ar° are required. Further- 
more, several elements in the ionized component - notably 
N, S - require small but important ionization corrections 
to infer accurate elemental abundances from the low-ion 
column densities. These corrections have important im- 
plications for our measurements of the metallicity of this 
damped system and our interpretations of the nucleosyn- 
thetic enrichment and dust depletion inferred from rela- 
tive chemical abundances. Before proceeding, we wish to 
emphasize that this system exhibits a number of charac- 
teristics (e.g. variations among the low-ion profiles) which 
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separate it from the majority of damped systems. There- 
fore, one must be careful not to generalize the ionization 
properties of this system to the entire sample of damped 
systems. At the same time, it is worth noting that the HI 
column density of this system is significantly greater than 
the damped Lya threshold of 2 x 10 20 cm~ 2 and would not 
have been expected to be partially ionized on the basis of 
its HI optical depth. 



Table 7 
METALLICITY TABLE 
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HM: One-Phase 
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KZ: Two-Phase, /(Fe+ 
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-1.36 
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These metallicity values are for the entire damped Lya system. The 
photoionization corrections, however, have only been applied to the 
gas at v ~ —26 km/s and —47 km/s. 

Consider first the impact of photoionization on the metal- 
licity measurements of this damped Lya system. We fo- 
cus on the metallicity of the neutral hydrogen gas because 
estimating the metallicity of all hydrogen gas (ionized + 
neutral) would require assumptions on the distribution of 
HI as a function of velocity. In addition, the metallicity of 
the neutral component is more readily compared against 
theoretical models. In Table 7, we report a series of metal- 
licity measurements for S, Si, and Fe assuming a variety of 
photoionization scenarios for the partially ionized gas 8 : (i) 
no corrections; (ii) a single-phase cloud with the HM spec- 
trum; (iii,iv) a two-phase scenario using the HM spectrum 
with two different assumptions on the level of ionization in 
the ionized component; (v-vii) single and two-phase sce- 
narios adopting the Kurucz spectrum with T = 40000 K. 
In the two-phase scenarios, the results are modestly sen- 
sitive to the AT(HI) value assumed for the ionized compo- 
nent. For this discussion, we assume iV(HI) = 10 19 cm~ 2 . 

The results presented in Table 7 indicate systematic off- 
sets in the metallicity of this damped system when one 
adopts ionization corrections. In general, however, the 
photoionization corrections imply a small, statistically in- 
significant decrease in the [X/H] values. The < 0.1 dex 
corrections are on the same order as the 0.1 dex uncer- 
tainty in the iV(HI) value. Therefore, even in the un- 
likely event that all of the damped Lya systems require 
ionization corrections which are this large, it will have 
a relatively minor impact on the chemical evolution his- 
tory currently resolved by the damped Lya systems (e.g. 
Prochaska & Wolfe 2000). The one possible exception, 
however, is the chemical enrichment history derived from 
Zn + whose photoionization balance is very poorly deter- 
mined (Howk & Sembach 1999). 

8 Note that the effects of photoionization on the metallicity of this 
system are tempered by the fact that only 50% of the gas requires 
ionization corrections. 
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Fig. 10. — The relative abundance [X/Fc] of the three main veloc- 
ity components in this damped system for the Haardt-Madau input 
spectrum. We consider (a) no ionization corrections; (b) a single- 
phase scenario; (c) a two-phase scenario with /jv(Fc+)= —0.3; and 
(d) a two-phase scenario with /jv(Fc+)= —0.1 dex. The x's mark 
the v fs —47 km/s component, the A's mark the v ~ —26 km/s 
component, and the d's mark the v ~ km/s component. 
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Fig. 11. — Same as Figure 10 but with the Kurucz input spectrum. 



In contrast to the metallicity measurements where fac- 
tors of 0.1 to 0.2 dex are generally unimportant, pho- 
toionization corrections can significantly influence our in- 
terpretations of the relative elemental abundances of the 
damped Lya systems. A correction to Si + /Fe + of 0.2 dex, 
for example, would imply a very different dust-to-gas ratio 
and/or nucleosynthctic enrichment history (e.g. Prochaska 
& Wolfe 2002). This point is well illustrated in Figures 10 
and 11 which present the relative abundance [X/Fe] of the 
three main velocity components in this damped system 
for the Haardt-Madau and Kurucz input spectrum respec- 
tively. In each figure we consider (a) no ionization correc- 
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tions; (b) a single-phase scenario; (c) a two-phase scenario 
with /jy(Fe + ) = —0.3; and (d) a two-phase scenario with 
/jv(Fe+) = -0.1 dex. 

Examining Figures 10 and 11 we find that the ratios 
of S/Fe and Si/Fe vary by > 0.1 dcx among the various 
ionization models. Therefore, it is very difficult to assess 
the elemental abundances ratios of the ionized gas in this 
damped Lya system. We are concerned that any damped 
system which is ionized to a similar degree will prove just 
as challenging. These uncertainties aside, we believe this 
system exhibits variations in its relative abundances inde- 
pendent of photoionization. These may be due to differen- 
tial depletion and/or varying nucleosynthctic enrichment 
patterns. In either case, this is one of the few examples 
of a damped Lya system where the abundances ratios are 
not uniform within statistical error. 

The results of this paper highlight the importance of ex- 
amining photoionization diagnostics in a larger sample of 
damped Lya systems. Although Al ++ /Si + may provide 
a rough assessment of ionization (Vladilo et al. 2001), a 
more precise analysis can be performed if one acquires ob- 
servations of adjacent ions like Fe + , Fe ++ and N°, N + . 
Of course, the intermediate-ions require observations with 
coverage deep within the Lya forest. While it may be 
important to consider a full photoionization analysis of 
the kind presented here for each damped Lya system, we 
believe the following photoionization diagnostics (summa- 
rized in Table 8) can provide a 'first-look' analysis of the 
ionization state. In all cases, our conclusions are derived 
from the extensive set of CLOUDY photoionization mod- 
els that we have described throughout the paper. 

• The Ar°/(Si + ,S + ) ratio is an excellent diagnostic 
for identifying systems with low ionization 
factors, x w 0. Specifically, a generic feature of 
photoionization models is that an observed value of 
[Ar°/(Si+,S+)] > -0.2 dex requires x < 0.1. For 
the 7 systems with Ar/Si measurements or limits 
(Molaro et al. 2001; Prochaska & Wolfe 2002), 
they all suggest x < 0.1 except GB1759+7539. 

It is also important to emphasize that a low 
Ar°/(Si + ,S + ) ratio does not require x » given 
its large reionization cross-section (Sofia & Jenkins 
1998). In fact, we have argued that the gas at 
v w km/s in the GB1759+7539 damped system 
has x < 0.1 even though it exhibits a significantly 
sub-solar Ar°/Si + ratio. Therefore, this diagnostic 
is only valuable for ruling out photoionization 
effects in a given damped system. Finally, we 
should note that the Ar°/(Si + ,S + ) ratio could be 
affected by nucleosynthesis and/or dust depletion. 
On the other hand, Ar/S is solar in systems with 
a wide range of metallicity in the local universe 
(Henry & Worthey 1999) and the majority of 
damped systems exhibit [S/Si] w indicating 
differential depletion will not modify the Ar, S, or 
Si abundances. 

• Unlike the Ar°/(Si + ,S + ) ratio, a determination 
of the Fe ++ /Fc + ratio can establish whether a 
given damped system is predominantly neutral 
or ionized. We find x < 0.1 in systems where 
Fc ++ /Fc+ < — 1.6 dex, and similarly, a system with 



Fc++/Fc+ > -1 dex will have x > 0.5. A further 
advantage of this diagnostic is that the Fe ++ /Fe + 
ratio is independent of nucleosynthesis or dust 
depletion. In short, an accurate measurement of 
Fc ++ via the Fc III 1122 transition is an extremely 
valuable probe of the ionization state of any quasar 
absorption line system. 

• The N + /N° ratio places similar constraints on x 
as the Fe ++ /Fc + diagnostic. We find x < 0.1 in 
systems showing N + /N° < — 1 dex and x > 0.5 in 
systems where N+/N > -0.2 dex. 



Table 8 

PHOTOIONIZATION DIAGNOSTICS 



Ratio 


x < 0.1 


x > 0.5 


Comment 


Ar u /(Si+,S+) 


> -0.2 


N/A 


Sensitive to dust, nuclco 


Fe++/Fc+ 


< -1.6 


> -1 


Lya forest 


N+/N 


< -1 


> -0.2 


Lya forest 



In passing, we wish to comment on the implications of 
photoionization for our kinematic modeling of the damped 
Lya systems (e.g. Prochaska & Wolfe 1997). To assess 
the kinematic characteristics of an individual damped Lya 
system we choose a single, high S/N low- ion profile as rep- 
resentative of the neutral gas for the entire system. This 
approach was motivated by the fact that the low-ion pro- 
files track one another very closely (e.g. Prochaska & Wolfe 
1996). In the case of the damped Lya system studied here, 
we have found significant variations from low-ion to low- 
ion even prior to photoionization corrections. This varia- 
tion will have a very small impact on the measured velocity 
width (Av; c.f. Prochaska & Wolfe 1997) of the damped 
Lya system, but the variations could impact the various 
shape parameters implemented in our analysis. Although 
this causes us some concern, we emphasize that this sys- 
tem is unusual for exhibiting a significant variation among 
the low-ion profiles. In a future paper, we will demonstrate 
the excellent chemical and ionization uniformity observed 
in the large majority of damped Lya systems (Prochaska 
et al. 2002). 

We also wish to reflect on the impact of photoionization 
regarding the relationship between the observed C II* pro- 
file and star formation (see Wolfe et al. 2002 for a more 
complete discussion). As noted in the introduction, the 
CII* profile of the GB1759+75 damped system does not 
closely track the other low-ion profiles (Figure 1). Specifi- 
cally, the C II* profile exhibits much greater optical depth 
at » K km/s than v « —26 km/s. We believe the gas 
at v « km/s has properties similar to the Cold Neutral 
Medium of the Galactic ISM (e.g. Savage and Scmbach 
1996). In contrast, if the gas at v w —26 km/s has phys- 
ical properties similar to the Warm Neutral Medium or 
Warm Ionized Medium of the Galactic ISM, then we would 
predict a low C II* optical depth relative to the other low- 
ions, as observed. Therefore, we believe the components 
at v « and —26 km/s represent two distinct phases of 
gas. The net effect is that it is difficult to assess the star 
formation rate within this damped Lya system from its 
measured C II* column density. 
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Finally, wc wish to outline other avenues for investigat- 
ing the ionization state of the damped Lya systems. One 
might gain important insight into the properties of ion- 
ized gas associated with damped systems by performing a 
comprehensive analysis of the quasar absorption line sys- 
tems with iV(HI) « 10 19 - 5 cm" 2 (e.g. Prochaska & Wolfe 
1999; O'Meara et al. 2001). These 'sub-DLAs' may show 
a wide range of ionization states and would provide an ex- 
cellent test of the modeling implemented throughout this 
paper. Furthermore, they are far more numerous than the 
damped systems and one may design specific observational 
experiments to focus on their ionization states (e.g. pre- 
select systems where Si ++ , N + etc. readily observed). An- 
other more efficient approach would be to examine lower 
rcdshift systems where the Lya forest is a far less com- 
promising effect. The challenge here remains identifying 
enough systems toward very bright quasars; a challenge 
which will be alleviated in part with the introduction of 
COS on the Hubble Space Telescope. 
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